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ABSTRACT 


Many countries in the world have opened their power markets to allow the competition 
between the power producers. The affect of the load changes on the operating frequency 
and real power output of the generator have been studied in the present work. If the 
frequency goes out of the limits the generator may loose synchronism .The frequency of 
the system should be maintained with in acceptable limits. In an interconnected power 
system, change in load of one control area will also affect the frequency and generation 
of the generators in other control areas. The variation in generation cost because of these 
deviations in frequency and real power output of the generator have been studied for a 
single control area case. The excessive generation cost because of the load change in this 
control area can be charged to that particular user. The powerful Newton-Raphson 
method has been used for solving the load flow problem. The real power generation has 
been optimized using optimal power flow. Load frequency control has been used to 
calculate the frequency deviations. 
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Chapter 1 


INTRODUCTION 


1.1 General 

Small scale power generating technologies are gradually replacing conventional 
generating technologies in various applications, in the electric power system. These 
distributed technologies have many benefits, such as high fuel efficiency, short 
construction lead time, modular installation and low capital expense, which all contribute 
to their growing popularity [1]. In the last few years many countries in the world have 
opened their power markets, to allow competition between power producers and at the 
distribution side also. The industry restructuring process is moving the power sector in 
general away from the traditional vertical integration and cost-based regulation toward 
increased exposure to market forces. 

The considerations involved in the smooth integration of distributed generation into the 
distribution system range from long term siting questions to concerns over maintaining 
frequency stability and the desired voltage profile. Once location and mode of operation 
are decided, and the necessary protection equipment is installed, the small generators will 
be able to supply power to customers, whether by contracting directly with customers, a 
power marketer or the system operator. The modeling for these bulk interactions involves 
using well established static models such as load flow and optimal flow models. 

The load flow solution gives the nodal voltages and phase angles and hence the power 
injection at all the buses and power flows through interconnecting transmission lines. 
Load flow solution also gives the initial conditions of the system when the transient 
behavior of the system is to be studied [2]. The Newton-Raphson method is used for 
solving the load flow problem. 

The optimal power flow uses the load flow solution and gives the solution, which results 
in optimal cost of generation [3]. 
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Research is going on how the users who cause the frequency changes at the generator can 
be penalized. Much work has been done in the previous years on the reactive power 
pricing. This thesis work is started to attain some price model to price the frequency 
changes at the generator. The load in the power system goes on changing time to time. 
Following a load change the frequency deviations and the generation changes at every 
generator are to be calculated. Some kind of control should be there to attain normal 
frequencies. In the competitive market there can be number of independent producers 
who will take decisions so as to maximize the profits and to maintain the operating 
voltage and frequency at the normal value. 

Various power systems [4] have been taken for the study. Some of these are ill- 
conditioned systems as given in [4]. 

1.2 Thesis Organization 


The Ne-wton-Raphson Load Flow method has been explained to calculate the voltage and 
phase angle at all buses for the given load schedule without constraints. Then by taking 
constraints also the load flow is run. The optimal power flow to optimize the generation 
obtained through load flow solution has been explained in chapter 2. Various test systems 
have been taken for the study and load flow solutions for these systems with and without 
constraints have been shown in chapter 2. 

When the load at one bus is increased the frequency and real power generation will also 
change at the generator. All the generators in every system have been taken as a single 
control area. Using the load flow method the load at the generator for the new condition 
has been calculated from which the change in the frequency and real power output of the 
generator can be calculated. The price model that calculates the change in the generation 
cost because of the change in the frequency and generation is explained in chapter 3. The 
generation cost for every test system has been given in chapter 3 considering every 
system as an isolated control area. 
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Chapter 2 


LOAD FLOW AND OPTIMAL POWER FLOW 
SOLUTIONS FOR VARIOUS TEST SYSTEMS 


2.1 INTRODUCTION 

2.1.1 Load Flow Solution 


Load flow solution is a solution of the network under steady state condition 
subjected to certain inequality constraints under which the system operates [2]. These 
constraints can be in the form of load bus voltages, reactive power generation of the 
generators, the tap settings of a tap changing under load transformer etc. 

The load flow solution gives the bus voltages and phase angles and hence the 
power injection at all the buses and power flows through interconnecting transmission 
lines. Load flow solution is essential for designing a new power system and for planning 
extension of the existing one for increased load demand. These analyses require the 
calculation of numerous load flows under both normal and abnormal (outage of 
transmission line, or outage of some generating source) operating conditions. Load flow 
solution also gives the initial conditions of the system when the transient behavior of the . 
system is to be studied. Single phase representation is enough under balanced operating 
conditions. A load flow solution of the power system requires mainly the following steps: 

• Formulation of the network equations. 

• Suitable mathematical technique for solution of the equations. 

Under steady state condition the network equations will be in the form of simple 
algebraic equations. The load and hence generation are continually changing in a real 
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power system, but for solving load flow it is assumed that loads and hence generation are 
fixed at a particular value over a suitable period of time. E.g. half an hour or so. 

In a power system each bus or node is associated with four quantities, real and 
reactive powers, bus voltages magnitude and its phase angle. In a load flow solution two 
out of four quantities are specified and the remaining two are required to be obtained 
through the solution of the equations. The buses are classified depending upon the 
quantities specified into the following three categories 

Load bus: At this bus the real and reactive components of power are specified. It is 
desired to find out the voltage magnitude V and phase angle 5 through the load flow 
solution. Voltage at load bus can be allowed to vary within the permissible value e.g. 5%. 

Generator bus or voltage controlled bus: Here the voltage magnitude corresponding to 
the generation voltage and real power Pq corresponding to its ratings are specified. It is 
required to find out the reactive power generation Qo and the phase angle 6 of the bus 
voltage. 

Slack, swing or reference bus: Here the voltage magnitude V and phase angle 6 are 
specified. This will take care of the additional power generation required and 
Transmission losses. It is required to find of real and reactive power generations (Pq, Qq) 
at this bus. 

Load flow solution can be achieved with iterative methods. There are many kinds 
of iterative methods out of which Newton-Raphson method is superior. The second order 
method is also used in the work for solving the load flow problem. 
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2.1.2 Optimal Power Flow (OFF) 


In the load flow problem as explained above, two variables are specified at each 
bus and the solution is then obtained for the remaining variables. The specified variables 
are real and reactive powers at PQ buses, real powers and voltages at PV buses, and 
voltage and angle at the slack bus The additional variables to be specified for load flow 
solution are the tap settings of regulating transformers. If the specified variables are 
allowed to vary in a region constrained by practical considerations (upper and lower 
limits of real and reactive generations, bus voltage limits and range of transformer 
settings), these result in infinite number of load flow solutions each pertaining to one set 
of values of specified variables. The best choice in some sense of the values of specified 
variables leads to the best load flow solution. Operating economy is naturally 
predominant in determining the best choice. 

The main economic factor in power system operation is the cost of generating real 

power. 
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2.2 NEWTON-RAPHSON LOAD FLOW (NRLF) METHOD 


2.2.1 Calculation of Jacobian 


For a N-bus power system there will be n equations for real power injection P, and 
n-equations for reactive power injection . 

Pd. = t Ik 'P. K, )] 

7=1 

e, = Qa. - Co = I [y. V, y, -s,-e , )] ; = 1,2,3, n 

/ = ! 

....( 2 . 1 ) 

The number of equations to be solved depends upon the specifications we have. If 
the total number of buses is n and number of generator buses is m then the number 
equations to be solved will be number of known P,’s and number of known Q,’s. In the 
above conditions number of known P,’s are n-1 and the number of known Q,’s are (n-m), 
therefore the total number of simultaneous equations will be 2*n-m-l, and number of 
unknown quantities are also 2*n-m-l. Unknowns to be find out are power angles (6) at all 
the buses except slack (i.e. n-1) and bus voltages (V) at load bus (i.e. n-m). The following 
mej;hod known as Newton-Raphson method is used for solving the unknown quantities. 
The problem formulation is as follows: 


(dP dP\ 


^ AP^ 


dS 

dV 








Ud 


\d5 dVj 
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Ai^ = P^ {specified ) - 
^Q:=Qi{specified)-Q^ 


Real power terms will be calculated for all the buses except slack bus and reactive 
power terms will be calculated for all load buses. In the above equation 


fAS\ 

Ud 


is the mismatch vector 

is the correction vector 


and 


J = 


dS dV 
dS dVj 


is the Jacobian matrix 


(2.4) 


The elements of the Jacobian matrix can be calculated using the following equations 


dP, 


dS, 

dP, 


dS, 

dP 


= V,V,Y„sm{6,-5,-e,^,) 


dV. 


= S rj,oos(s,-s,-e„)+v,G„ 


/ /-I 


dV. 


V,Y„cos{d,-d^-e,) 


P -V^ G 

j. I r j 


^Q, ^ n rr2 

acj, 

35 . 


v,vj,iCos(s,-s,-e„) 


V^Y,sin{s.~S,-e„Yv,B., 
|^ = r,y„sin(5,-<5 ,-^,) 


.(2.5) 
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Procedure for this iterative method is for the given system first the Y-bus matrix has to be 
formed. 

Y=G+jB 

Where 

Y is admittance matrix 
G is real part of Y-bus matrix 
B is imaginary part of Y-bus matrix 

The resistance and reactance of each line have been given for the given system with 
which the admittance matrix can be calculated. 


2.2.2 Iterative Algorithm for N-R Method 

1 . With voltage and angle (usually J = 0 ) at slack bus fixed, assume voltage 
magnitude and power angles at PQ buses and S at all PV buses. Generally flat voltage 
start will be used. 

2. Compute AP, for all buses except slack bus and Ajg, for all PQ buses 

using Eq. (2.3). If all the values are less than the prescribed tolerance, stop the iterations. 

3. If the convergence criterion is not satisfied, evaluate elements of the 
Jacobian using Eq. (2.5) 

4. Solve the Eq. (2 2) for correction vector. 

5. Update voltage angles and magnitudes by adding the corresponding 
changes to the previous values and return to step 2. 

2.2.3 Adjust ments in NRLF Solution 

In the present thesis work some constraints have been applied on voltages at the 
PQ buses, reactive power at the PV buses so as to keep the system healthy. 
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2.2.3.1 Constraints on reactive power: 


Since Q, at generator buses is not given, it is calculated in each iteration at all generator 
buses and it is been checked for the condition 

a™.„<a<amax -(2.6) 

If the above mentioned condition satisfied then the i* bus will remain as generator 
bus and there will be no change in the procedure. 

If the condition violates then Q, will be set at the limit values it can be either lower limit 
or upper limit depending upon the violation and i* bus will be treated as load bus from 
the next iteration. Because of this change of bus type one additional equation 
corresponding to AQ, will be added to NRLF scheme. 

2.2.3.2 Constraints on bus voltages; 

After every iteration voltages and angles are updated with the obtained correction 
matrix. Voltage has to satisfy the conditioned given below 

V <V <V r? 7i 

If it satisfies the condition then there will be no change in the process, if it violates then 
voltage has to be fixed at one of its lower or upper limits depending upon the violation 
and the bus type is been changed to generator type. Because of this one equation 
corresponding to AQ, will be deleted from NRLF scheme, so the size of Jacobian matrix 
will be decreased. 
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2.2.4 Flow Chart for NRLF method 
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Fig 2.1 Flow chart for load flow solution using Newton-Raphson method 







2.3 OPTIMAL POWER FLO W 


2.3.1 Generator Operating Cost 


The major component of generator operating cost is the fuel input/hour, while 
maintenance contributes only to a small extent. The input-output curve of a unit can be 
expressed in million kilocalories per hour or directly in terms of rupees per hour versus 
output in megawatts. For every generator, is the minimum loading limit below 

which it is uneconomical to operate the unit and is the maximum output limit. 

The anal 5 l;ical operating cost can be written as C, (P^., ) Rs/hour at output , where the 

suffix i stands for the unit number. It generally suffices to fit a second degree polynomial 
[3], i.e., 

C, Rs/hour (2.8) 


The slope of the cost curve, i.e., — ^is called the incremental fuel cost (IC), and 


clP. 


Gt 


is expressed in units of rupees per megawatt hour (Rs/MWh). If the cost curve is 
approximated as in Eq. (2.8) then the expicssion for IC will be 

{lC),=a,P,:,+b, (2.9) 

i.e., a linear relationship. 


2.3.2 Optimal Load Flow Sol nf ion 


The solution technique is based on load flow solution by NR method, a first order 
gradient adjustment algorithm for minimizing the objective function and use of penalty 
functions to account for inequality constramls on dependent variables. 
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2.3.2.1 OPF Without Inequality Constraints 


The objective function to be minimized is the operating cost 

C = EC,(-Po,) 


.( 2 . 10 ) 


subject to the load flow equations 


p.-ZK\ 


7=1 

n 


2: + I 

7=1 


V. 


V. 


Y, 

cos ( 

0, + 


)=0 " 


sin ( 

Ptj + 

-Sj. 

)=0 



J 


> for each PQ bus 


and 


P, -l]|F,||l>;,||i;|cos(^y +S, -Sj)=0 foreachPVbus (2.11) 

/=i 


It is to be noted that at the z'* bus 


P = P -P 

•* 1 ^Gi ^Di 


Qi ~Qgi Qdi 

Where Po, and Qp, are load demands at bus z . 

Equations (2.1 1) can be expressed in vector form 


.( 2 . 12 ) 


f(x,y)= 




7=1 

n 


7=1 

7=1 


r For each PQ bus 


For each PV bus 


.(2.13) 
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Where the vector of dependent variables is 


X = 



For each PQ bus 
For each PV bus 


and the vector of independent variables is 


(2 14) 


r. 



For slack bus 
For each PQ bus 

For each PV bus 


u 


y = 


p 


(2.15) 


In the above formulation, the objective function should include the slack bus power. 

The vector of independent variables y can be partitioned into two parts - a vector u of 
control variables, which are to be varied to achieve optimum value of the objective 
function and a vector p of fixed or disturbance or uncontrollable parameters. Control 
parameters may be voltage magnitudes on PV buses, at buses with controllable power 
etc. 

The optimization problem can now be restarted as 


min C(x,u) 


(2.16) 


subject to equality constraints 

f(x,u,p) = 0 


(2.17) 


To solve the optimization problem, defining the Lagrangian function as 
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(2.18) 


£(x,u,p) = C(x,u)+X ^ f (x,u,p) 

where 1 is the vector of Lagrangian multipliers of the same dimension as f(x,u,p). 

The necessary conditions to minimize the unconstrained lagrangian function are 

d£ dc raff, ^ 

dx dx Lax_ 

a£ dC [dfT, ^ 

au au L5u_ 

|^ = f(x,u,p) = 0 (2.21) 

Where f, x, u, p and 1 are vectors. 

The Eq. (2.21) is obviously the same as equality constraints. The expressions for 

aC df . 

and "Tyas needed in equations (2.19) are rather involved. — is nothing but 

ax au dx 

Jacobian matrix. 

Equations (2.19), (2.10) and (2.21) are non-linear algebraic equations and can 
only be solved iteratively. A simple yet efficient iteration scheme, that can be employed, 
is the steepest descent method (also called gradient method). The basic technique is to 
adjust the control vector u, so as to move from one feasible solution point in the direction 
of steepest descent (negative gradient) to a new feasible solution point with a lower value 
of objective function. By repeating these moves in the direction of the negative gradient, 
the minimum will finally be reached. 

23.2.2 Com f utational Procedure 

The computational procedure is as follows. 

1 . Make an initial guess for u, the control variables. 

2. Find a feasible load flow solution by NR iterative method. The method 
successfully improves the solution x as follows 


( 2 . 20 ) 


(2.19) 
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where Ax is obtained by solving the set of linear equations 

df 


or 


dx 


(x('»,y) Ax = -f(x(''),y) 


Ax = -(jM)’'f(xMy) 

the end solution of step 2 are a feasible solution of x and the Jacobian matrix. 
3. Solve the Eq. (2. 1 9) for 




- 1-1 


ax 


.( 2 . 22 ) 


4. Insert X from Eq. (2.22) into Eq. (2.20), and compute the gradient 

.nr 


V£ = 


ac 


au 


K 

du 


.(2.23) 


af 


It may be noted that for computing the gradient, the Jacobian J = — is already known 

ax 

from the load flow solution. 

5. If V£ equals zero within prescribed tolerance, the minimum has been 
reached. Otherwise, 


6. Find a new set of control variables 

(2-24) 

where 

Au = -aV£ (2.25) 


Here Au is a step in the negative direction of the gradient. The step size is adjusted by 
the positive scalar a . 

Selection of a is the critical part of the algorithm. Too small a value of 
a guarantees the convergence but slows down the rate of convergence; too high a value 
causes oscillations around the minimum. 
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23.2.3 Inequality Constraints on Control Variables 


The control variables are always constrained 


u ■< u < u 

“min — “ — “max 


e.g. 

Pr <Pr ^Pr ...( 2 . 26 ) 

G/,min ~ G/ G/,max 

If the correction Aw, in Eq. (2.25) causes w, to exceed one of the limits, w, is set equal to 
the corresponding limit, i.e.. 


r 

U 

/,max 


U 


i,new 




^/,mm 


if u At/ > u 
if u + Aw < U 

J I, old ^ ^/,mm 

otherwise 


(2.27) 

After a control variable reaches any of the limits, its component in the gradient 
should continue to be computed in later iterations, as the variable may come within limits 
at some later stage. 

The necessary conditions for minimization of £ under constraint Eq. (2.26) are: 


T- = 0 if <W, <W,,.nax 

OM, 


du, 


<0 


f 


u, =u 


I, max 


5m, 


>0 


f u, 


■ u 


/,mm 


(2.28) 


Therefore, in step 5 of the computational algorithm, the gradient vector has to satisfy the 
optimality condition Eq. (2.28). 
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Once the generation has been decided for the given load schedule the load flow program 
will be run again after some period as the demand goes on varying. In a power system 
when the load disturbance occurs the operating frequency and real power output of the 
generator will change. To determine the change in the load at each generator for the new 
load schedule again the load flow will be run. These load changes will be the input 
vector to the load frequency control problem to determine the deviations in the frequency 
and real power output. 
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2.4 RESULTS 


43 bus [4], 13 bus [4], IEEE 14 bus and 11 bus [4] test systems have been taken for the 
study. The data is given in the appendix A. The Newton Raphson Load Flow method has 
been used to solve the load flow problem. The load flow problems have been solved 
without constraints and with constraints. The optimal power flow has been run to 
optimize the generation cost. The results of the NRLF method and OFF are shown in the 
following sections. 


2.4.1 Results for 43 bus test system 

The load flow problem for 43 bus test system [4] has been solved without putting any 
constraints. The general base case solution is shown in the Table 2.1. 


/ 
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Table 2. 1 Load flow solution for 43 bus test system [4] using NRLF method 


without any constraints 


Voltage 
(P-«) 

Phase Angle 
(deg.) 

1.136000 

0.000000 

1.038553 

-11.196495 

1.015844 

-15.235573 

0.984997 

-13.712234 

1.035090 

-11 313513 

1.037332 

-11.568977 

1.040850 

-15.315901 

1.009215 

-13.401521 

1.007399 

-13.422747 

1 008570 

-13 487796 

1.008570 

-13.487796 

1.043339 

-15.553347 

0 986969 

-13.673939 

1.009973 

-14 641201 

1.034337 

-11.451182 

0.972859 

-17.680355 

1.008570 

-13.487796 

1.006835 

-15.975134 

1.032887 

-16.428598 

1.045197 

-14.405982 

0.998752 

-13.240570 

1.029823 

-16.484206 

1.006628 

-13.304210 

1.022492 

-17.118298 

0.995748 

-13.431604 

1.010820 

-21.370186 

1.020050 

-15.432568 

1.031474 

-11.617084 

1.003496 

-13.178171 

0.919032 

-11.951977 

1.000000 

0.692896 

1.000000 

-6.065698 

1.000000 

-16.209620 

0.845448 

-21.630065 

1.000000 

-15.702249 

0.920797 

-21.756725 

0.975098 

-4.542215 

0.908746 

-21.630065 

1.026546 

-11.738752 

0.913340 

-12.175063 

1.005302 

-14.698162 

0.884260 

-19.544946 

1.024592 

-11.763712 
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The load flow problem has been solved using the first order NRLF method by 
considering constraints on bus voltages. The minimum and maximum limits have been 
taken as 0.9 p.u and 1.1 p.u. The solution is given in the Table 2.2. 


Table 2.2 Load flow solution using NRLF method for 43 bus test system [4] with the 

constraints on bus voltages 


Voltage 

(p.u) 

Phase Angle 
(deg.) 

1.136000 

0.000000 

1.100000 

-10.553851 

1.100000 

-14.290340 

1.052288 

-12.973546 

1.100000 

-10.711639 

1.100000 

-10.904547 

1.100000 

-14.274729 

1.072088 

-12.558681 

1.070422 

-12.579715 

1.071666 

-12.644764 

1.071666 

-12.644764 

1.100000 

-14.467805" 

1.053497 

-12.906119 

1.083747 

-13.703580 

1.100000 

-10.845451 

1.041819 

-16.320639 

1.071666 

-12.644764 

1.079117 

-14.917460 

1.100000 

-15.236020 

1.100000 

-13.485208 

1.061412 

-12.368501 

1.100000 

-15.374980 

1.069123 

-12.444268 

1.100000 

-15.759599 

1.060151 

-12.609425 

1.100000 

-19.577074 

1.100000 

-14.352148 

1.100000 

-11.058515 

1.065570 

-12.300383 

0.926179 

-11.152750 

1.000000 

0.747890 

1.000000 

-5.312052 

1.000000 

-15.323243 
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0.847790 

-20.728674 

1.000000 

-14.817287 

0.923281 

-20.854644 

0.998606 

-4.363649 

0.911263 

-20.728674- 

1.095469 

r -11.168957 

0.920555 

-11.372754 

1.075776 

-13.760233 

0.893243 

-18.610630 

1.093655 

-11.191561 


After getting the load flow solution with constraints on bus voltages the optimal power 
flow is run for the system to optimize the generation cost. The real power output of each 
generator is limited to a maximum value of (1.1*Pgi) p.u and to a minimum value of 
(0.9*Pg,) p.u. The solution giving the bus voltages and phase angles at each bus after 
optimal power flow considering the limits on the real power output of the generator is 
shown in Table 2.3. 


Table 2.3 Voltage magnitudes and Phase Angles at each bus after Optimal Power Flow 

for 43 bus test system [4] with constraints 


Voltage (p.u) 

Phase Angle (deg.) 

1.136000 

0.000000 

1.100000 

-9.787865 

1.100000 

-12.938932 

1.048144 

-11.749151 

1.080971 

-9.665371 

1.099894 

-10.102128 

1.100000 

-12.919196 

1.064307 

-11.284862 

1.062623 

-11.305931 

1.063859 

-11.370980 

1.063859 

-11.370980 

1.100000 

-13.308905 

1.049021 

-11.676514 

1.083086 

-12.910483 

1.099966 

-10.064411 
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1.033322 

- 15.105379 

1.063859 

- 11.370980 

1.074265 

- 13.705566 

1.100000 

- 14.013776 

1.100000 

- 12.440043 

1.055842 

- 11.114266 

1.099511 

- 14.126218 

1.062189 

- 11.181230 

1.100000 

- 14.182029 

1.055060 

- 11.368999 

1.100000 

- 17.908845 

1.100000 

- 12.944973 

1.099240 

- 10.262341 

1.059518 

- 11.048646 

0.925502 

- 9.940663 

1.000022 

2.063675 

1.000036 

- 4.097875 

1.000001 

- 14.172171 

0.847557 

- 19.559920 

1.000003 

- 13.657221 

0.923035 

- 19.685959 

0.996339 

- 3.054983 

0.911013 

- 19.559920 

1.094849 

- 10.372876 

0.919872 

- 10.160957 

1.075143 

- 12.967147 

0.892393 

- 17.411172 

1.093034 

- 10.395500 
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2.4.2 Results forl3 bus test system 


The load flow problem for 13 bus test system [4] has been solved without any constraints. 
The general base case solution is shown in the Table 2.4. 


Table 2.4 Load flow solution for 13 bus test system [4] using NRLF method 

without any constraints 


Voltage (p.u) 

Phase Angle (deg.) 

1.000000 

0.000000 

1.063717 

1.347255 

1.179639 

1.945324 

1.160323 

1.984012 

1.000000 

2.110602 

1.037000 

9.211764 

1.063684 

8.365741 

1.100000 

7.426077 

0.943000 

13.840178 

1.100000 

7.826597 

0.957134 

11.835802 

1.036184 

7.363695 

, 0.969457 

4.326548 
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The load flow problem has been solved using the first order NRLF method by puttiing 
constraints on bus voltages. The minimum and maximum limits have been taken as 0.9 
p.u and 1 . 1 p.u. The solution is given in the Table 2.5. 


Table 2.5 Load flow solution using NRLF method for 13 bus test system [4] with the 

constraints on bus voltages 


Voltage (p.u) 

Phase Angle (deg.) 

1.000000 

0.000000 

1.063578 

1.387019 

1.100000 

2.247327 

1.100000 

2.247327 

1.000000 

2.248149 

1.037000 

9.473620 

1.063746 

8.679546 

1.100000 

7.807725 

0.943000 

14.221826 

1.100000 

8.208245 

0.957134 

12.217450 

1.036184 

7.745343 

0.969457 

4.708196 
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After getting the load flow solution with constraints on bus voltages the optimal power 
flow is run for the system to optimize the generation cost. The real power output of each 
generator is limited to a maximum value of (l.l*Po,) p.u and to a minimum value of 
(0.9*Pgi) P-u. The solution giving the bus voltages and phase angles at each bus after 
optimal power flow considering the limits on the real power output of the generator is 
shown in Table 2.6. 


Table 2.6 Voltage magnitudes and Phase Angles at each bus after Optimal Power Flow 

for 1 3 bus test system [4] with constraints 


Voltage (p.u) 

Phase Angle (deg.) 

1.000000 

0.000000 

1.064659 

1.239276 

1.100000 

1.987836 

1.099962 

1.987937 

1.000038 

1.988586 

1.040593 

8.466641 

1.062673 

7.757496 

1.092275 

6.959252 

0.967112 

12.347970 

1.100000 

7.084377 

0.989823 

10.584633 

1.030016 

6.784220 

0.928995 

4.103847 
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2.4.3 Results for IEEE 14 bus test system 


The load flow problem for IEEE 14 bus test system has been solved without any 
constraints. The general base case solution is shown m the Table 2.7. 


Table 2.7 Load flow solution for 14 bus test system using NRLF method 

without any constraints 


Voltage (p.u) 

Phase Angle (deg.) 

1.060000 

0.000000 

1 045000 

-4.517426 

1 010000 

-11.962797 

1.021593 

-9.367419 

1.027582 

-7.874867 

1.070000 

-11.655461 

1.062376 

-12.008487 

1.090000 

-12.008487 

1.055376 

-13.383275 

1.050294 

-13.365262 

1.056302 

-12.648834 

1.055290 

-12.580481 

1.050081 

-12.733563 

1.034994 

-14.102923 
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The load flow problem has been solved using the first order NRLF method by puttiing 
constraints on bus voltages. The minimum and maximum limits have been taken as 0.9 
p.u and 1 .1 p.u. The solution is given in the Table 2.8. 


Table 2.8 Load flow solution using method for IEEE 14 bus test system with the 

constraints on bus voltages 


Toltage(i).«) 

Phase Angle (deg.) 

1.0600DO 

0.000000 

1 019702 

-4.209865 

0.980821 

-12.077397 

0.988329 

-9.248864 

0.995728 

-7.651070 

1.00094-3 

-11.613080 

1.018032 

-12.217016 

1.047298 

-12.217016 

1 .004776 

-13.786931 

0.995999 

-13.738793 

0.994520 

-12.850663 

0.98661 3 

-12.671585 

0.982425 

-12.880401 

0.975829 

-14.525790 
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After getting the load flow solution for the problem with constraints on bus voltages the 
optimal power flow is run for the system to optimize the generation cost. The real power 
output of each generator is limited to a maximum value of (1.1 *Pg.) p.u and to a 
minimum value of (0.9*Pgi) p.u. I he solution giving the bus voltages and phase angles at 
each bus after optimal power flow considering the limits on the real power output of the 
generator is shown in Table 2.9. 


Table 2.9 Voltage magnitudes and Phase Angles at each bus after Optimal Power Flow 
for IEEE 14 bus test system with constraints 


Voltage (p.u) 

Phase Angle 
(deg.) 

1.060000 

0.000000 

1.044759 

-4.517426 

1.009980 

-11.962797 

1.021593 

-9.367419 

1.027582 

-7.874867 

1.070003 

-11.655461 

1.062376 

-12.008487 



1.090001 

-12.008487 

1.055376 

-13.383275 

1.050294 

-13.365262 

1.056302 

-12.648834 

1.055290 

-12.580481 

1.050081 

-12.733563 

1.034994 

-14.102923 


29 




2.4.4 Results for 11 bus test system 


The load flow problem for 1 1 bus test system [4] has been solved without any constraints. 
The method has diverged for a tolerance level of 0.0001. It has converged when the 
tolerance is taken as 0.008. But the values obtained are not practical. Practically 1 1 bus 
test system diverging when there are no constraints on voltages. 

The load flow problem has been solved using the first order NRLF method by 
considering constraints on bus voltages. The minimum and maxitnnm limits have been 
taken as 0.9 p.u and 1 . 1 p.u. The solution is given in the Table 2. 1 1 . 


Table 2.1 1 Load flow solution using NRLF method for 1 1 bus test system [4] with the 

constraints on bus voltages 


Voltage (p.u) 

Phase Angle (deg.) 

1.024000 

0.000000 

1.066178 

-2.903824 

1.056163 

-4.662305 

1.041775 

-3.852663 

1.044332 

-5.537236 

1.061867 

-3.934891 

0.900000 

-24.872891 

0.900000 

-23.957400 

0.900000 

-23.878411 

0.900000 

-36.085608 

0.900000 

-37.681380 
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Chanter 3 


GENERATION PRICING IN SINGLE AREA 

POWER SYSTEM 


3.1 INTRODUCTION 


Power being a complex quantity, it has been observed that any change in real 
power output of the generators directly affect the system frequency and tie-line loadings 
while any change in reactive power affect only the system voltage. This property helps in 
dividing the control problem of a power system into two separate channels: the MW- 
frequency control channel and the M VAR- voltage control channel. The control of the 
real power output of the generators in response to change in system frequency and the tie- 
line power flow so as to maintain the scheduled system frequency and the tie-line 
interchange with in permitted limit is termed as Automatic Generation Control (AGC) or 
Load Frequency Control (LFC) [6]. On the other hand the control of reactive power 
generation in the system in order to maintain constant voltage is known as excitation 
Control. 

The main objective of the AGC for an inter connected power system can be stated 
as follows 

1 . Matching generation to load 

2. Regulating system frequency error to zero 

3. Distributing generation amongst areas so that inter-area tie-line flows 
match a prescribed schedule. 

4. Distributing generation with in each area such that the operating cost of 
area is minimized. 
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Automatic Generation Control therefore can be subdivided into two separate 
control problems. One is, the traditional Load Frequency Control problem, which meets 
the first three objectives. Second is, Economic Load dispatch problem, which takes care 
of the last objective. 

Generally objective (1) is achieved by the system governor. However, system 
governor alone is found inadequate to take care of objectives (2) and (3). Therefore, 
supplementary control is added to the system governor utilizing PID controller such that 
the deviations in frequency and tie-line loading from the prescheduled values, following a 
sudden load change in any area, are reduced to zero in the time span of less than a 
minute. The controller design should be such that not only the objectives (2) and (3) 
should be met but also the transient oscillations in frequency are kept to a minimum . 

In a competitive market there will be independent generators in a large number. Every 
producer will take his own decisions while supplying power to the corresponding 
consumers. He will fix the cost per unit generation depending upon the conditions in 
which the generators are being operated and the changes in the demand. In every aspect 
he will try to increase his profit and at the same time he will keep in mind the competition 
from the other producers. Industry restructuring, and particularly the deregulation of 
generation, is opening the power sector to market forces. 

To coordinate distributed generator actions in the short-term operations and 
control of the spot energy and ancillary services markets, a price model is proposed. The 
method used is developed to accurately capture the cost associated with local deviations 
from the scheduled power and energy. A mathematical framework for price model, 
designed to coordinate distributed generators as they participate in both the short run 
energy market and the ancillary services market. This mathematical framework will take 
care of both frequency change as well as load change. 
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3.2 LOADJBEC UENCY CONTROL 


3.2.1 Model of Governor 


The speed governing system can be modeled as in Fig 3.1 


AP,{s) 



Fig. 3.1 Model of Speed Governor 


Where 

R=speed regulation of the governor 
K =gain of speed governor 

=time constant of speed governor 

A/),, is the commanded increase in power. This signal sets into motion a sequence 
of events. When the load increases the frequency goes down. As a result turbine 
generator speed decreases. Then this signal comes into picture to increase the steam input 
thus increasing speed, which in turn increases, the frequency. 

AT is the output signal of the governor, i.e., input signal to the turbine. 

3.2.2 Model of Turbine 


The model requires a relation between changes in power output of the steam 
turbine to changes in its steam valve opening. A non-reheat turbine with a single gain 
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fflctor K~i snci s single time constant is considered. Thus the trunsfer function of the 
turbine is 



_£l_ 

1 + 7] i' 


(3.1) 


The block diagram representation is shown in Fig 3.2 


Ay(.0 


^ 




l + T,s 

^ 



Fig. 3.2 Model of Turbine 

3.2.3 Model of Generator and Load 

The increment in power input to the generator-load system is 

AP,-AP, 

where AFfi = APi , incremental turbine power output and APy is the load increment. 

As the frequency changes, the motor load changes being sensitive to speed, the 
rate of change of load with respect to frequency, i.e., dPp/8f can be regarded as nearly 
constant for small changes in frequency A/ and can be expressed as 
(dPi,/^)Af = BAf 

where the constant B can be determined empirically. B is positive for a motor load. 
Writing the power balance equation, 

AP„-AP„-^|-(A/)+SA/' (3.2) 

AP„(pu)-APa{pu)==^~iAf)+B(jm)Af (3.3) 

Taking Laplace transform, 
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(3.4) 




D 2// 

B + -~- s 

r 


AFis)=[AP^.{s)-AP,{s)] 


K 


ps 


+ s j 


.(3.5) 


where 


T. 


2H 


1 


=Powcr system time constant 


~~B sain 


The block diagram representation is shown in Fig 3.3 


APa 



3.2.4 Complete Block Diagram for Single Area Power System 

The complete block diagram representation of an isolated power system 
comprising turbine, generator, governor and load is easily obtained by combining the 
block diagrams of individual components. The complete block diagram with feedback 
loop is shown in fig (3.4). 
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3.4 Block diagram representation of Load Frequency Control 



Two important incremental inputs to the load frequency control system are - , 

the change in speed changer setting; and , the change in load demand. If a simple 
situation in which the speed changer has a fixed setting (i.e., AP(.=0) and the load 
demand changes. This is known as free governor operation. For such an operation the 
steady state change in system frequency for a sudden load change by an amount APyj is 




is steady state frequency 


(3.7) 
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3.2.6 Dynamic Res )onse 


1 0 obtain the dynamic response giving the change in frequency as function of the 
time for a step change in load, the Laplace inverse of Eq. (3.6) is to be calculated. 

The droop in the frequency from no load to full load should be as small as 
possible as much change in frequency cannot be tolerated. In fact the steady change in the 
frequency should be zero. While steady state frequency can be brought back to the 
scheduled value by adjusting speed changer setting, the system could undergo intolerable 
dynamic frequency changes with changes in load. It leads to the natural suggestion that 
the speed changer setting be adjusted automatically by monitoring the frequency changes. 
For this purpose, a signal from A/’ is fed through a PID controller to the speed changer. 
The above modification makes the steady state error to fall to zero and settling time to 
reduce. 

Therefore, 

A/). (.v) ={-K^- Kjs - K,s)b AF{s) (3.8) 

where b is feed back gain. 

A', is gain of proportional controller 
A'j is gain of integral controller 
A", i,s gain of differential controller 

In this case the change in frequency will be fed back to the PID controller the output of 
the output of the controller is A?( .(,y). 
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I'ig. 3.5 Load Frequency Control with PID controller 

The signal A/».(.v) generated by the PID controller must be opposite sign to AF{s) which 
accounts for negative sign in the block for integral controller. Now 

.,{i + t,4\ + r,sXi + + -J + 

~ £f.vAF(s) = 0 (3.9) 


3.3 GENERATION PRICING 

With in the distributed system load is distributed throughout the system and 
generators are located at specific buses. To simulate the dynamic behavior of the syst 
disturbances are specified as load fluctuations of small magnitude to allow th 
small-signal, liner models. A system model is defined by specifying the dis 
system topology, the location and size of loads and the location, size and yp 
generators. The inputs to the models are the system disturbances, represented as the input 
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vector to the system of state equations. The output from the simulation is the dynamie 
behavior of all the state variables, with frequency and real power output typically being 
of greater interest than the others. Monitoring the frequency at each bus can assess the 
frequency stability. 

The goals in developing models for analyzing frequency behavior are to represent 
the dynamics ol' distributed generators in response to system disturbances, and to propose 
and analyze the effectiveness of different control strategies designed to ensure system 
stability. 

'I’he modeling effort is based on building deeoupled, linearized state space models 
for each distributed generator and coupling them through a distribution system model. 
The models that include a synchronous generator all use a form of the swing equation [1] 
as the generator stale equation: 

J5+D8 = P„-P^ ( 3 . 10 ) 

Where 1\. s , the electrical power output. 

is the mechanical power from the turbine. 

'fhe turbines in the system are assumed as steam turbines to continue the 
following analysis. In developing the model of the generator, the objective is to represent 
each generator with a small number of state variables. 

The small signal dynamic for each generator [1] is 

MWf. *(e, -£>)W(. +/1 -Pc, 

T,P, ^-P,+k,a 

7;a = -w,,~ra (3 11^ 

Where 

M is the inertia constant 

e, is a coefficient representing the turbine self-regulation defined as dP, / dwg 
D is the damping coefficient 
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7^^ is the time constant representing the delay between the control valves and the 
turbine nozzles 

ki is a proportionality factor representing the control valve position relative to the 
turbine output variation 

is the time constant of the valve-servomotor-turbine gate system 

a is the gate position 

r is the permanent speed droop of the turbine 
After getting the deviations in the frequency and generation in each control area 
following the load changes in every area, the generations cost in every control area are 
calculated as follows. 

The price model presented is for decoupled real power/frequency d 3 mamics. The 
reason for this is two fold. First, with this emphasis, the modeling effort mirrors the 
pattern to date for developing a spot price or responsive price system, which usually 
focuses on pricing real power, since that is the major commodity of the industry. The 
second reason is that the use of distributed generators for voltage support is reasonably 
well accepted by the power industry. In contrast, the frequency dynamics of distribution 
systems with distributed generation units, and the possibility of these units participating 
in the supply of ancillary services such as frequency stability and spinning reserve, are 
relatively new issues. Small signal, linearized models are used for analyzing these 
markets. 

I’he development of the price model begins here by expressing the cost of power 
generation in terms of the state variables in the generator equations. Cost can be 
incorporated into the state space generator models by writing an output equation to 
capture the variable costs associated with generating power from any given technology. 
Each state space model identifies the set of elements that together can reproduce the basic 
machine performance. The cost equation [1] then will become now, 

cost = C =C^WQ+CpPi + C^a + C{PQ) (4.12) 

where ^ is given by the Eq. (2.8) 
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The coefficients in the above equation represent the marginal cost associated with each 
piece of equipment or process represented by the specified state variable. 

With the addition of the output cost equation, the model for each generator can be 
expressed as 


Mw,;={e^-D)wc+P,-Pa 

T„ P; = -P, +k, a 

T^a = -WG-ra + w,^ (3.13) 

cosf = C =C^WQ + CpP,-vC,a + C{Pc) 

The above model can be represented as 


C = h(x,PG) 


(3.14) 


The generators and the system will respond to the price signal at specific 
intervals, indicating that the price signal is best modeled in discrete time. The first step in 
developing this discrete time model is to assume the primary dynamics have settled. 
Therefore the liq. (3.14) will now become 


C = /z(x,P,) 


(3.15) 


Solving these equations for cost results in a discrete time cost equation of the form 


C{k)=y^ WQ{k)+C{PQ{k))+r2W^^f (3.16) 


Where k represents the discrete time index for the price control loop, and and y^ 
are constant c.xprcssions of the generator parameters and cost coefficients. For the steam 
generator, these coefficients are of the from 


Yf 

Yi 


Cpk, c; 




V 


J 


(3.17) 
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T4 RESULTS 


All the generators have been assumed to be in a single control area. It is assumed 
that all the generators are having similar characteristics such that these generators can be 
taken a single ecjuivalcnt generator oi the same type. The marginal costs are taken as 
shown in the 1 able A. 7. As the real system data is not available the marginal costs 
associated with each piece of equipment have been assumed. 

Change in load at any bus will ultimately increase load on generators, which will 
cause frequency deviation at the generator. Current price model has been used to 
calculate the generation cost considering both the deviation in frequency and the change 
in real power output. Various test systems have been taken for the study. Data for these 
systems is shown in appendix A. For a step load change at one of the buses the load flow 
has been run to find out the actual load change that will occur at the generator, this 
change has been taken as AP^ (.v) and frequency deviation and generation change have 
been calculated using the method explained in section 4.2. Using these quantities 
generation pricing has been done as explained in section 4.3. Results are shown below. 

3.4.1 Generation Pricing for 43 bus test system 

43 bus test system shown in Fig. (A.l ) is taken for the study. The data is given in 
the appendix A. 'I'he generator model parameters [1] are shown in Table. A. 8. 

The load at the bus 3 has been increased by 0.05 p.u. The load flow program is 
again run and the actual load burden on the generators have been calculated which has 
been taken as A/^,. 

By using Controller this frequency deviation can be brought back to zero, here 
they have not been used. The main reason for not using the controllers is the users that 
are drawing excessive load from the scheduled power themselves will reduce the usage of 
the power as the price is increasing because of the deviations in real power output and 
frequency. As the demand is more than the supply the frequency settled down at a lower 
value. The final frequency, generation and price are tabulated below before disturbance, 
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m pu 


After getting the deviations in frequency and real power generation the actual 
frequency and generation have been calculated. The base frequency is taken as 50 Hz. 
I'he frequency plot is shown in Fig. 3.8 and the total generation is shown in Fig. 3.9. 
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Fig. 3.8 Total generation in 43 bus power system after a load change of 0.05 p.u at bus 3 
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Fig. 3.10 Generation cost in 43 bus power system after a load change of 0.05 p.u at bus 3 

Table ^.1 Comparison of the frequency, real power output and generation cost before and 

after disturbances for 43 bus test system 


Quantity 

Before disturbance 

After disturbance 

Peak value 

I’requcncy (I Iz) 

50 

49.7591 

49.6073 

Cieneration (p.u.) 

11.5047 

11.585 

11.6167 


2.1143*10^ 

2.1746*10^ 

2.1913*10^ 





3.4.2 Generation Pricing in 13 bus test system 


1 3 bus test system [4] shown in Fig. (A.2) is taken for the study. The data is given 
in the appendix A. The generator model parameters [1] are shown in Table. (A.8). 

The load at the bus 12 has been increased by 0.02 p.u. The load flow program is 
again run and the actual load burden on the generators have been calculated which has 
been taken as AP^ . 

By using Controller this frequency deviation can be brought back to zero, but here 
they have not been used. The main reason for not using the controllers is the users that 
are drawing excessive load from the scheduled power themselves will reduce the usage of 
the power as the price is increasing because of the deviations in real power output and 
frequency. As the demand is more than the supply the frequency settled down at a lower 
value. The final frequency, generation and price are tabulated below before disturbance, 
after disturbance and their peak values in table 4.2. The deviation in frequency is shown 
in Fig. 3.11. The change in generation is shown in Fig. 3.12. 






After getting the deviations in frequency and real power generation the actual 
frequency and generation have been calculated. The base frequency is taken as 50 Hz. 
The frequency plot is shown in Fig. 3.13 and the total generation is shown in Fig. 3.14. 



I'ig. 3.13 Total generation in 13 bus power system after a load change of 0.02 p.u at bus 

12 
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Th. generation cost is calculated using the price model after getting the deviations in the 
frequency and the real power output. The variation in the cost is shown in Fig.3, 15, 
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Fig. 3.15 Generation cost in 1 3 bus power system after a load change of 0.02 p.u at bus 

12 


Table 3.2 Comparison of the frequency, real power output and generation cost before and 

after disturbances for 13 bus test system 


Quantity 

Before disturbance 

After disturbance 

(steady state) 

Peak value 

Frequency (Hz) 

50 

49.9443 

49.9092 

Generation (p.u.) 

1.8 

1.8584 

1.8657 

Price (Rs/hr.) 

1.5543*10^ 

1.5557*10^ 

1.5575*10^ 
















3-4.3 Generation Pricing in 14 bus test system 


14 bus test system shown in Fig. (A.3 ) is taken for the study. The data is given in 
the appendix A. The generator model parameters [1] are shown in Table. (A.8). 

The load at the bus 9 has been inereased by 0.05 p.u. The load flow program is 
again run and the actual load burden on the generators have been calculated which has 
been taken as A/^, . 

By using Controller this frequency deviation can be brought back to zero, but here 
they have not been used. The main reason for not using the controllers is the users that 
are drawing excessive load from the scheduled power themselves will reduce the usage of 
the power as the price is increasing because of the deviations in real power output and 
frequency. As the demand is more than the supply the frequency settled down at a lower 
value. The final frequency, generation and price are tabulated below before disturbance, 
after disturbance and their peak values in table 4.3. The deviation in frequency is shown 
in Fig. 3.16. The change in generation is shown in Fig. 3.17. 
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After getting the deviations in frequency and real power generation the actual 
frequency and generation have been calculated. The base frequency is taken as 50 Hz. 
1 he frequency plot is shown in Fig. 3.18 and the total generation is shown in Fig. 3.19. 



f'ig. 3.18 Total generation in 14 bus power system after a load change of 0.05 p.u at bus 9 
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The generation cost is calculated using the price model after getting the deviations in the 
frequency and the real power output. The variation in the cost is shown in Fig.3.20. 
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Fig. 3.20 Generation cost in 14 bus power system after a load change of 0.05 p.u at bus 9 


Table 3-3 Comparison of the frequency, real power output and generation cost before and 

after disturbances for 14 bus test system 


Quantity 

Before disturbance 

After disturbance 

Peak value 

Frequency (Hz) 

50 

49.8695 

49.7873 

Generation (p.u.) 

2.74 

1 

2.7944 

2.8116 

Price (Rs/hr.) 

2.1708*10'^ 

2.4545*10"^ 

2.4842*10"^ 
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3.4.4 Generation Pricing in 11 bus test system 


1 1 bus test system [4] shown in Fig. (A.4) is taken for the study. The data is given 
in the appendix A. The generator model parameters [1] are shown in Table. (A.8). 

The load at the bus 5 has been increased by 0.05 p.u. The load flow program is 
again run and the actual load burden on the generators have been calculated which has 

been taken as APjr, . 

By using Controller this frequency deviation can be brought back to zero, but here 
they have not been used. The main reason for not using the controllers is the users that 
are drawing excessive load from the scheduled power themselves will reduce the usage of 
the power as the price is increasing because of the deviations in real power output and 
frequency. As the demand is more than the supply the frequency settled down at a lower 
value. The final frequency, generation and price are tabulated below before disturbance, 
after disturbance and their peak values in table 4.4. The deviation in frequency is shown 
in Fig. 3.21. The change in generation is shown in Fig. 3 22. 


61 









pu generatioi 



Fig. 3.22 Change in Generation of 1 1 bus power system after a load change of 0.05 p.u at 





After getting the deviations in frequency and real power generation the actual 
Ircquency and generation have been calculated. The base frequency is taken as 50 Hz. 
The frequency plot is shown in Fig. 3.23 and the total generation is shown in Fig. 3.24. 


Actual Generation 



Fig. 3.23 I’otal generation in 1 1 bus power system after a load change of 0.05 p.u at bus 5 
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The generation cost is calculated using the price model after getting the deviations in the 
frequency and the real power output. The variation in the cost is shown in Fig.3.15. 


Price plot 



Fig. 3.25 Generation cost in 1 1 bus power system after a load change of 0.05 p.u at bus 5 


'fable ^.2^Comparison of the frequency, real power output and generation cost before and 

after disturbances for 13 bus test system 


Quantity 

Before disturbance 

After disturbance 

Peak value 

Frequency (Hz) 

50 

49.849 

49.75 

Generation (p.u.) 

0.7708 

0.8215 

0.8414 

Price (Rs/hr.) 

3158 

3237 

3261 
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Chapter 4 


CONCLUSIONS AND SCOPE FOR FURTHER 

WORK 


4,1 Conclusions 

The load flow problems for various test systems have been solved using NRLF method 
with and without constraints. The general base case solution for 1 1 bus system without 
constraints is not a better solution as the voltage magnitudes are very fer from the 
operating voltage. Where as if the limits on the bus voltages are considered it has 
converged to a better solution. For 13 bus system also if the constraints are not taken the 
voltages are higher in magmtude. When the constraints are considered the voltage 
magnitudes have been limited with in the prescribed limits. For the other systems there is 
no such difference in the solution even when the constraints are considered. 

In Optimal Power flow the objective fimction has been taken as the minimi/ing 
the generator operating cost. For every system the optimal power flow has been run to 
reschedule the generation, which is obtained from the load flow solution to reduce the 
generation cost. 

All the systems have been considered as single isolated control areas. When the 
load at any bus changes it will cause deviation in system frequency as well as real power 
generation. This study has been done on four test systems. In every system the deviations 
in frequency and real power generation have been calculated following a load change at 
any bus. The generation cost has been taken as a fimction of both frequency and real 
power generation. When the load increases the generation cost will increase and vice- 
versa as the frequency decreases and generation increases. The PID controllers have not 
been used in the thesis work. The reason is when any user draws more than scheduled 
power the increase in generation cost will be charged to that particular user who will then 
shed the load automatically. Then the frequency will settle back to the normal value. 
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4.2 Scope for Further work 




• The work can be extended to penalize the users those are causing frequency 
changes by drawing unscheduled power in such a way this extra charge which 
they have to pay will make the load to come down and the frequency to go 
upto the normal value. 

• The operating cost can be minimized by taking the constraints on operating 
frequency also. Where as in the present thesis work only the load flow 
constraints have been taken. 

• In this thesis work every system has been taken as a single control area. One 
can take a system such that it can be divided into a number of control areas 
and then the price model can be extended for calculating generation cost in 
every area following the load changes. 
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Appendix A 


Data For 43 Bus Test System 



116 MW 
53 MVR 


3« MW 
30 MVR 
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The 43 bus system is shown in Fig. A.l. The system data is taken from ref. [4, 7], 
The relevant data are provided in following tables. Table A.l gives Y Bus matrix 
elements of 43 bus system. Table A.2 gives the operating condition of the system. 


Table A.l Y Bus matrix elements of 43 bus system 


From 
Node (i) 

To Node 

(j) 

S13 

Bx, 

1 

1 

0 

- 30.609 

1 

2 

0 

30.609 

2 

2 

481.288 

“ 1545.194 

2 

5 

- 277.195 

873.583 

2 

6 

- 34.368 

108.124 

2 

15 

“ 169.726 

534.322 

3 

3 

0 

“ 5.714 

3 

4 

0 

6.015 

4 

4 

61.331 

“ 69.160 

4 

13 

“ 61.331 

62.874 

5 

5 

277.195 

“ 916.892 

5 

7 

0 

21.277 

5 

8 

0 

20.513 

6 

6 

34.368 

“ 118.699 

6 

12 

0 

10.638 

7 

7 

0 

-20 

8 

8 

452.840 

“ 482.861 

8 

9 

“ 288.938 

295.777 

8 

23 

“ 163.902 

167.191 

9 

9 

300.983 

- 317.044 

9 

10 

“ 12.045 

12.342 

9 

16 

0.0 

8.796 

10 

10 

12.045 

“ 20.885 

10 

11 

0.0 

2.857 

10 

17 

0.0 

5.714 

11 

11 

0.0 

“ 2.857 

12 

12 

0.0 

- 10.0 

13 

13 

92.381 

“ 100.709 

13 

18 

0.0 

6.015 

13 

25 

- 31.05 

31.640 

14 

14 

0.0 

“ 15.015 

14 

43 

0,0 

15.4 

15 

15 

340.398 

“ 916.783 

15 

19 

0.0 

8.649 

15 

20 

0.0 

15.791 

15 

28 

“ 170.673 

357.003 

16 

16 

0.0 

“ 8.576 

17 

17 

0.0 

“ 5.714 

18 

18 

0 

“ 5.714 

19 

19 

164.292 

“ 280.783 
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19 

22 

- 164.292 

272.805 

20 

20 

^ 0 

- 15.002 

21 

21 

104.312 

“ 143.609 

21 

24 

0 

9.267 

21 

29 

- 104.312 

133.623 

22 

22 

164.292 

- 282.281 

22 

26 

0 

9.023 

23 

23 

321.579 

- 328.81 

23 

29 

- 157.677 

161.76 

24 

24 

0 

- 8.572 

25 

25 

87.15 

- 106.814 

25 

27 

0 

9.023 

25 

29 

- 56.1 

65.824 

26 

26 

0 

- 8.572 

27 

27 

0 

- 8.572 

28 

28 

373.447 

- 612.837 

28 

39 

- 202.775 

256.136 

29 

29 

318.089 

- 372.311 

29 

30 

0 

3.766 

29 

37 

0 

7.895 

30 

30 

125.789 

- 524.464 

30 

32 

0 

30.769 

30 

38 

0 

4.131 

30 

40 

- 125.789 

485.547 

31 

31 

0 

- 13.038 

31 

37 

0 

13.038 

32 

32 

0 

- 30.769 

33 

33 

0 

- 3.320 

33 

38 

0 

3.32 

34 

34 

0 

- 7.365 

34 

38 

0 

6.852 

35 

35 

0 

- 6.18 

35 

38 

0 

6.18 

36 

36 

0 

- 2.703 

36 

38 

0 

2.703 

37 1 

37 

0 

- 21.348 

38 

38 

0 

- 22.398 

39 

39 

512.581 

- 663.260 

39 

41 

0 

15,015 

39 

43 

- 309.806 

392.255 

40 

40 

125.789 

- 508.837 

40 

42 

0 

21.622 

41 

41 

0 

- 15.015 

42 

42 

0 

-20 

43 

43 

309.806 

- 408.029 
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Table A.2 Operating condition of 43 bus system 


Bus 

No. 

Voltage 

magnitude 

(g-^) 

Phase 

Angle 

(deg.) 

Net Real 
Power (p.u) 

Net Reactive 
Power (p.u) 

1 

1.136 

0 

0 

0 

2 

1 

0 

0 

0 

3 

1 

0 

“ 0.16 

“ 0.12 

4 

1 

0 

0 

0 

5 

1 

0 

“ 0.53 

- 0.4 

6 

1 

0 

0 

0 

7 

1 

0 

“ 1.6 

- 1.2 

8 

1 

0 

0 

0 

9 

1 

0 

0 

0 

10 

1 

0 

0 

0 

11 

1 

0 

0 

0 

12 

1 

0 

- 0.8 

“ 0.6 

13 

1 

0 

0 

0 

14 

1 

0 

“ 0.8 

- 0.6 

15 

1 

0 

0 

0 

16 

1 

0 

“ 0.64 

“ 0.48 

17 

1 

0 

0 

0 

18 

1 

0 

“ 0.24 

“ 0.18 

19 

1 

0 

0 

0 

20 

1 

0 

- 0.88 

“ 0.66 

21 

1 

0 

0 

0 

22 

1 

0 

0 

0 

23 

1 

0 

0 

0 

24 

1 

0 

“ 0.64 

“ 0.48 

25 

1 

0 

0 

0 

26 

1 

0 

- 0.8 

“ 0.6 

27 

1 

0 

- 0.32 

“ 0.24 

28 

1 

0 

0 

0 

29 

1 

0 

0 

0 

30 

1 

0 

0 

0 

31 

1 

0 

1.16 

0.52 

32 

1 

0 

2.9 

2.57 

33 

1 

0 

0.285 

0.3 

34 

1 

0 

0 

0 

35 

1 

0 

0.58 

0.56 

36 

1 

0 

- 0.005 

0.030 

37 

1 

0 

0 

0 

38 

1 

0 

- 1.44 

“ 1.02 

39 

1 

0 

0 

0 

40 

1 

0 

0 

0 

41 

1 


“ 0.8 

“ 0.3 

42 

1 

0 

“ 2.24 

“ 1.68 

43 

1 

0 

0 

0 
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Table A.3. Bus data for 13 bus system 


Bus 

No. 

Bus 

Voltage 

(pu) 

Phase 
Angle 
(deg. ) 

Generation 

Load 

Real 

Power 

(pu) 

Reactive 

Power 

(pu) 

Real 

(pu) 

Reactive 

(pu) 

1 

1 . 0* 

0.0 


- 

1.65 

0.56 

2 

1.0 

0.0 

0 

0 

0 

0 

3 

1.0 

0.0 

0 

0 

0 

0 

4 

1.0 

0.0 

0 

0 

0 

0 

5 

1.0 

0.0 

0 

0 

0 

0 

6 

1.037* 

0.0 

0.5 

- 

0.05 

0.03 

7 

1.0 63 

0.0 

0 

0 

0 

0 

8 

1.1* 

0.0 

0 

- 

0 

0 

9 

0.943* 

0.0 

0.5 

- 

0 

0 

10 

1 . 1* 

0.0 

0 

- 

0 

0 

11 

1.0 

0.0 

0 

0 

0.05 

0.05 

12 

1.0 

0.0 

0 

0 

0.05 

0.032 

13 

1.0 

0.0 

0 

0 

0 

0 


* Input data 


Table A.4. Line data for 13 bus ill-conditioned system 


Branch 

From 

To 

Resistance 

Reactance 

Susceptance 

Number 

Node 

Node 

(pu) 

(pu) 

(pu) 

1 

1 

2 

0.0040 

0.085 

0 

2 

1 

3 

0.0040 

0.0947 

0 

3 

5 

4 

0.0040 

0.0947 

0 

4 

4 

3 

0.0074 

0.1430 

0.436 

5 

6 

2 

0.0481 

0.4590 

0.246 

6 

6 

7 

0.0090 

0.1080 

0.016 

7 

8 

3 

0.0121 

0.233 

0.712 

8 

7 

8 

0 ^ 

0.15 

0 

■ ■ 9’ "■ 

9 

10 

0.0105 

0.2020 

0.620 

10 

10 

11 

0 

-0.15 

0 

11 

11 

12 

0.0086 

0.1665 

0.508 

12 

12 

13 

0.0075 

0.1465 

0.448 

13 

13 

8 

0 

-0.15 

0 


Base=1000 MV A 
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Table A. 5. Transformer data for 13 bus system 


Branch number 

From node 

To node 

Tap setting 

1 

1 

2 

+ 5% 

2 

2 

3 

+ 10% 

3 

5 

4 

+ 10% 


Table A. 6 Assumed Values for the Cost Coefficients 


Coefficient 

Value 

a 

1.4 

b 

0.8 

d 

10 


Table A.7 Marginal costs associated with various equipments 


Coefficient 

Value 

Cp 

175 


200 


175 


Table A. 8 Generator Model Parameters 


Parameter 

Value 

M 

1.26 

D 

2.0 

et 

0.15 

Tu 

0.2 

kt 

0.95 

Tg 

0.25 

R 

2.4 

. 
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DATA FOR 14 BUS SYSTEM 


The 14 bus system is shown in Fig A.3, data is taken from ref. [9] and buses are 
renumbered. The relevant data are provided in following tables. Table A.9 gives 
Generator Data, Table A. 10 gives Generator Bus voltage, Table A. 11 gives Transformer 
Data, Table A. 12 gives Load Bus data and Table A. 13 gives Line Data. 



Fig A.3 IEEE 14 bus test system 
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Table A.9: Generator Data for 14 bus system 


Generator 

No 

Real Power Generation 
Limit 

Reactb 

General 

/e Power 
-ion Limit 

Maximum 

(MW) 

Minimum 

(MW) 

Maximum 

(MW) 

Minimum 

(MW) 

1 

200.00 

050 0 

100.0 

-45.0 

2 

100.0 

020.0 

050..0 

-40.0 

3 

- 

- 

040.0 

00.0 

6 

150.0 

020.0 

024.0 

-06.0 

OO 

- 

- 

024.0 

-06.0 


Table A. 10: Generator Bus voltage Data for 14 bus system 


Bus 

No. 

Scheduled Real 
Power Generation 
Pg(MW) 

Specifieed Voltage 
Magnitude 

V spec(P-U) 

Load 

Real 

(MW) 

Reactive 

(MVAR) 

1 

- 

1.060 

00.00 

00.00 

2 

40.0 

1.045 

21.70 

12.70 

1 

3 


1.010 

94.20 

19.00 

6 

20.0 

1.070 

11.2 

7.5 

8 

- 

1.090 

00.00 

00.00 
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Table A.l 1 : Transformer Data for 14 bus system 


Line 

From 

To 

Series Impedence 

Tap Setting 

No 

Bus 

Bus 

Resistance (O.u) 

Reactance (p.u) 

8 

4 

7 

0.0 

0.2091 

0.978 


4 

9 

0.0 

0.5561 

0.969 


5 

6 

0.0 

0.2502 

0.962 


Table A. 12 Load Bus Data for 14 bus system 


Bus no 

Load 

External Shunt 
Susceptance 
(p.u) 

Real 

Reactive 

4 

47.8 

4.0 

0.0 

5 

7.6 

1.6 

0.0 

7 

0.0 

0.0 

0.0 

9 

29.5 

16.6 

0.19 

10 

9.0 

5.8 

0.0 

11 

3.5 

1.8 

0.0 

12 

6.1 

1.6 

0.0 

13 

13.5 

5.8 

0.0 


14.9 

5.0 

0.0 


Base MVA=100 
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Table A. 13: Line Data for 14 bus system 


Line. 

No. 

From 

No 

To 

Bus 

Series Impedance 

Shunt 

Susceptance 

(p.u) 

Resistance (p.u) 

Reactance (p.u) 

1 

1 

2 

0.01938 

0.05917 

0.528 

2 

1 

5 

0 05403 

0 22304 

0.0492 

3 

2 

3 

0.04699 

0.19797 

0.0438 

4 

4 

4 

0.05811 

0.17632 

0.0374 

5 

2 

5 

0.05695 

0.17388 

0.0340 

6 

3 

4 

0.06701 

0.17103 

0.0346 

7 

4 

5 

0.01335 

0.04211 

0.0 

11 

6 

11 

0.09798 

0.19890 

0.0 

m 

6 

12 

0.12291 

0.25581 

0.0 

13 

6 

13 

0.06615 

0.13027 

0.0 

14 

7 

8 

0.0 

0.17615 

0.0 

15 

7 

9 

0.0 

0.11001 

0.0 

16 

9 

10 

0.03181 

0.08450 

0.0 

17 

9 

14 

0.12711 

0.27038 

0.0 

18 

10 

11 

0.08205 

0.19207 

0.0 

19 

12 

13 

0.22092 

0.19988 

0.0 

20 

13 

14 

0.17093 

0.34802 

0.0 
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DATA FOR 11 BUS SYSTEM 


The 1 1 bus system is shown in Fig. A.4. The system data is taken from ref [4,7], 
The relevant data are provided in following tables. Table A. 14 gives bus data and Table 
A. 15 gives Y matrix elements 




Fig A.4 1 1 bus ill-conditioned system 
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Table A. 14 Bus data for 1 1 bus system 


Bus 

Number 

Voltage (p.u.) 

Angle 

(radians) 

Real power 
injection (p.u.) 

Reactive power 
injection (p.u.) 

1 

1.024 

0 

0 

0 

2 

1 

0 

0 

0 

3 

1 

0 

■a.i28 

-0.062 

4 

1 

0 

0 

0 

5 

1 

0 

-OA65 

-0.080 

6 

1 

0 

■0.090 

-0.068 

7 

1 

0 

0 

0 

8 

1 

0 

0 

0 

9 

^ 1 

0 

-0.026 

-0.009 

10 

1 

0 

■■■■■■I 

0 

11 

1 

0 


-0.057 


Table A. 1 5 Y bus matrix for 1 1 bus system 


From Bus i 

To bus j 

Conductance G,j 
(p.u.) 

Snsceptance Bij 
(p.u.) 

1 

1 

0.0 

-14.939 

1 

2 

0.0 

14.148 

2 

2 

12.051 


2 

3 

0.0 

6.494 

2 

4 

-12.051 

13.197 

3 

3 

2.581 


3 

5 

-2.5 81 

3.789 

4 

4 

12.642 

-74.081 

4 

5 

0.0 

2.177 

4 

6 

0.0 

56.689 

4 

7 

-0.592 

0.786 

5 

5 

2.581 

-5.889 

6 

6 

0.0 

-55.556 

7 

7 

3.226 

-4.304 

7 

8 

-2.213 

2.959 

8 

8 

2.893 

-5.468 

8 

9 

-.138 


8 

10 

-0.851 

1.163 

9 

9 


-1.042 

10 

10 

1.346 

-6.11 

10 

11 

-0.374 

3.742 

11 

11 


-2.785 


Base MVA=100 
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